% 读取文件，设置保留原始列名
data = readtable('附件3.xls', 'VariableNamingRule', 'preserve');

% 提取英国、美国、加拿大的 GDP 数据
UK_GDP = table2array(data(:, '英国'));
US_GDP = table2array(data(:, '美国'));
Canada_GDP = table2array(data(:, '加拿大'));

% 将数据合并成一个矩阵，每列代表一个国家的 GDP 时间序列
gdp_matrix = [UK_GDP, US_GDP, Canada_GDP];

% 截取 1980 - 2011 年的数据
% 假设数据是按年份顺序排列的，且每年有 1 个数据点
% 这里计算 1980 - 2011 年的行数
start_row = 1;
end_row = 2011 - 1980 + 1; 
gdp_matrix = gdp_matrix(start_row:end_row, :);

% 手动指定 VAR 模型的滞后阶数，这里假设为 2，你可以根据经验调整
aicOrder = 2;

% 建立 VAR 模型，明确变量个数为 3
model = varm(aicOrder, 3); 

% 估计 VAR 模型
results = estimate(model, gdp_matrix);

% 查看 VAR 模型的估计结果
disp(results);